Level statistics in a two-dimensional disc with diffusive boundary scattering 
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o\ : 

We calculate the energy level statistics in a two-dimensional disc with diffusive boundary scatter- 
ing by the means of the recently proposed ballistic nonlinear a- model [ffl. 
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1. The purpose of this paper is to use the recently proposed nonlinear a- model |1[ for ballistics in disordered 
conductors with long mean free path and find out how it works. In a recent paper || Andreev and Altshuler (AA) 
suggested a general method of calculation of level statistics in a disordered system beyond the limits of the random 
i— i . matrix theory. Their calculations were performed for a diffusive disordered system with the mean free path I shorter 
than the system size R. For this case and for energy difference u> exceeding the mean level spacing A (oj ^> A) their 
method was based on using the non-linear <7-model ||] and accounting for a perturbative contribution from the vicinity 
of several stationary points of the action M. AA have also conjectured a general form which the level statistics obeys. 
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In this paper we follow the same general strategy, addressing this problem for a quantum particle in a two-dimensional 

S disc with no scattering in the bulk and strong boundary scattering [Q. Such a problem is a natural target for the 
recently proposed field theory for quantum ballistics jj],™ . Consideration is restricted solely to unitary symmetry. In 
all relevant parts our results coincide with those independently obtained by Blanter, Mirlin and Muzykantskii [||. 

2. In order to apply the general approach of AA, we begin with replacing the non-linear diffusive super-matrix 
cr-model |?| by its ballistic generalization |lj . The partition function of this field theory is determined as a functional 
integral over a supermatrix g(n, r) = U~ 1 AU on the energy shell E = p 2 /2m, (n = p/p) in the phase space: 
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| where (..) = J ...dO/2w, dO = dO n ,dO' = dO a < and the scattering probability W(n, n') in the bulk is connected 
i with the mean free time r and transport mean free time r^ r as 
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In this paper we will consider a clean disc with no scattering in the bulk and strong boundary scattering. Therefore, 
both r and r^ r will be taken as being infinitely large. Integration in Eq (||) is not defined unless boundary conditions are 
imposed on super-matrix g at the inner boundary of the sample. Since super-matrix g has a meaning of a distribution 
function of electrons, the boundary condition it obeys is similar to that, which is applied to the distribution function in 
Q , classical kinetics R] || . General boundary condition for matrix-functions (see || for an example) is pretty complicated. 
Fortunately, for the purposes of this paper (calculation of the spectral correlation function with the precision to the 
first non-vanishing term beyond Random matrix theory) the problem could be significantly simplified, because most 
of the properties of the energy levels could be determined by the the values of g-matrix close to the special point 
<7(rn) = A. If U = 1 — w/2 + w 2 /8 + .., then the free energy F could be expressed through matrices w, what gives in 
- ' the quadratic approximation: 

Fo{w} = —-£- f drdO n sti[w 2 i(L- iu)wi2], (4) 

where indices 1 (2) relate to "retarded" ("advanced") degrees of freedom, and L denotes operator of the kinetic 
equation. Since the free energy in Eq (Q) is quadratic, it results in classical linear equation. The boundary condition 
which should be imposed upon 1012(1", n), is now a direct analog of the condition, imposed upon the distribution 
function in classical kinetics. Extremely strong boundary scattering is popularly modelled by the diffusive boundary 
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condition 0, @), which assumes that the distribution function for outgoing particle does not depend on angular 
variable n and is coupled to that for incoming particle by flux conservation. If J\f is an outward normal to the 
sample's boundary, then the diffusive boundary condition reads as 

w 12 {nM<0)= [ — nAAwi 2 (n'). (5) 

3. According to AA, the level statistics is determined by the determinant of a linear operator L from Eq (Q). The 
eigenvalue condition is 

nwfVw = Aw, (6) 

subject to boundary condition 

2w < = — w < / d4>cos(j)= / dcf> cos (jnu > (((>), (7) 

Jtt/2 J-tt/2 

where w < and «;> are the values of "distribution function" w(n) at the disc boundary at nAT > and nAf < 
respectively and w < does not depend on its argument. 

The left hand side of Eq (^J) consists of a derivative dw/dl along the trajectory of a particle inside the disc (see Fig 
1). Its solution has the form of a simple exponential 



w(l) = w(0) exp 



(8) 



Solution (^) should be substituted into the boundary condition (R). It is also convenient to express the direction 
of momentum cos <f> of incident electron at point 9 in of the disc boundary Eq (]?]) through that coordinate on the 
boundary 9' , where this electron was diffusively scattered from (cos 4> = sin[(0' — 9)/2\) . This all leads to the eigenvalue 
equation in the form: 
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The expansion of w < (9) in the Fourrier series w < (9) — ^2w m e" ne transforms the condition (|^) into 

i r 

fm(fJ>m,k) = 0, / m (/i) = 1 - - / exp [2imu + ^sinu] sinudu — 0, (10) 
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where \x m .k — 2R\ m .k/vp. One can see from Eq ( JTc| ) that one of the eigenvalues with m — vanishes ( say /io.o = 0). 
This corresponds to w independent of both n and r and it is not surprising that the relaxation rate of this eigen-mode 
vanishes. Substitution m — ► — m into Eq ( |l0| ) makes it clear that fj, m — /x_ m . The equation, complex conjugate to 
Eq ( |l0| ) shows that if \i is an eigenvalue, then /i* is an eigenvalue as well. None of the eigenvalues has a negative real 
part. A natural labeling || is k = 0, ±1,±2, ... for even m and k = ±1/2, ±3/2.. for odd m. For k = and even m 
the eigenvalues are real. The asymptote of the eigenvalues is 

lnfc / 1\ 

Mm.fe ~ — ^- + tt« I k + - I < m < k. (11) 

So, for < to -C k Im/i 3> Re[i and both don't depend on m. (Some the eigenvalues of the Liouvillcan are shown in 
Fig 2.) 

4. The purpose of this paper is to calculate the spectral correlation function 

R 2 {uj) = (7rAi? 2 ) 2 (i/(e)i/(e + u)) - 1, (12) 

where v(e) is the density of states and A = l/irR 2 v is the mean level spacing. The time of ballistic flight along 

diameter of the disc tf = 2R/vp introduces a natural scale for the frequencies. 

As has been shown by AA, the deviation of R2(u>) from the Wigner-Dyson expression 
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R2(s)=5(s)-—2-, s = - (13) 
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at frequencies w > A is well described by introducing the spectral determinant D(s) 

d(s) = n 
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which is closely connected with the spectral function S(uj), first introduced by Altshuler and Shklovskii |X(J for diffusive 
systems: 
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The spectral correlation function can be decoupled at ui 3> A into the sum p[ of a smooth part i? s 
oscillating part i? osc : 
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So, the calculation of the spectral determinant -D(s) is a key point of the whole AA programme, which we approach 
now. 

5. It is possible to write down an expression for the spectral determinant without an explicit computation of the 
eigenvalues. In order to do that, note that the function / m (/i), defined by Eq (10), is an entire function of its argument, 
which has only simple zeros at fj, = [im^ and f{\i)j{\i /(m)) vanishes as fi — > oo. Therefore, fmil^) can be represented 
as an infinite product (m ^ 0). 
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For m = the function /o(m) vanishes at /i = 0. So the same theorem could be applied to the function — 4/ (A t )/ 7r A t - 
Multiplying f m (fJ.) and f m (— (J,), taking the product over all m, and analytically continuing to /i = ±i£ = ±icutf, we 
arrive, finally, at the expression for the spectral determinant 
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where it is taken into account that 
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Since \ m ,k = fj. m ,k/tf, the spectral determinant -D(s) consists of two dimensionless parameters utf and At/. One of 
this parameters is always small (Atf -C 1), while the second one wtj could be either larger or smaller than unity. 
These two limiting cases constitute the limits of high and small frequencies respectively. 

6. At low frequencies < 1 the spectral determinant -D(£) can be simlified and the asymptotes of both the 
smooth and the oscillatory parts of the spectral correlation functions coincide, as was first discovered by Kravtsov 
and Mirlin Q. This gives the following expression for the spectral correlation function 
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Using the low frequency asymptote of Eq (^]), we can present the spectral function 

S(u>) = 5 m (w) 
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in the form of a contour integral 



Sm(co) = ^-<f> 
27TI Jr 
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where contour C encloses all zeros of of the function f m (z). As u> — > we obtain from Eq ( p2[ ) the following expression 
for the coefficient B in Eq (|2l| ) 
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7. In order to find the asymptote of the spectral determinant D in the high frequency limit, consider a product 
P(i£) 
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Its logarithm is presented by the sum 



ln_P(i£) = In 1 / dusmuexp(2imu + i£sinit) 
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At this tage, it is convenient to use the indentity 
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which replaces the sum over m by the sum over n and integral over x. For large values of £ the integral over u in 
Eq(^||) is small and the logarithm should be expanded up to the second order in this integral (linear term vanishes). 
After calculation the sum over n the expression could be simlificd to the following form: 
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To evaluate the spectral determinant, we need to find the product P(i£)P(— i£) Using the steepest descent method, 
we arrive at £ 3> 1 at following asymptote for the spectral determinant: 



This gives the smooth part of the spectral correlation function in the form: 
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Eq (|T^) gives the oscillatory part of the spectral function equal to 
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8. In conclusion, we found that the application of the ballistic non-linear a- model [[TJ to the study of level statistics 
for electrons in a clean disc with strong boundary scattering enables us to solve this problem beyond the limits of the 
random matrix theory. 

A clean disc with diffusive scattering on its boundaries, unlike other chaotic systems, has an upper limit for the 
time of flight at t = tf = 2R/vf- Therefore, if a Fourier transform of a time dependnet form-factor is calculated, 
it oscillates as a frequency w with period 2ir/tf. As it was shown in section 7, the smooth part of the spectral 
corerelation function i?2(w) at high frequencies uitf 3> 1 is proportional to square of the relevant form- factor. This 
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leads to oscillations of R sm (uj) with twice as shorter period ir/tf. In our understanding, such oscillations would not 
appear in a general case. 

Another striking result is exhibited in Eq (|3(]) : the amplitude of the oscillatory part of spectral correlation function 
is small as (Atf) 2 , but does not decay with lo, unlike one obtained by AA for a diffusive system. This could be 
understood if to recall that our disc is clean inside and, therefore, at short times t <C tf certain correlation between 
electron wave function remains. This correlation is small and proportional to (pfR)~ 2 ~ (Ai/) 2 , but decays with the 
energy difference to much slower. If our disc has a bulk disorder with the meen free time r 3> tf Jl2[ , R osc oc exp[— lot]. 
To similar result leads variation of the Fermi velocity with energy: R osc oc exp[— uj/Ep]. The smooth part of spectral 
correlation function (see Eq ( |29| ) ) also exhibits weak dependence on energy difference to. In our understanding, these 
results are of a general nature. 
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